摘 要:本数据集基于2013‒2018年NPP/VIIRS夜间灯光数据,进行影像合成、不稳定光源去除以及逐年校正等处理,消除短暂性光源和背景噪声的影响。将已处理的夜间灯光数据与人口统计数据进行相关性分析建立回归模型,按照误差分类进行人口空间化建模,以消除山地地形、人口数量少等因素对模拟精度的影响,最终生成海南岛人口空间分布数据集。该数据集以.tif格式存储,空间分辨率为500 m,由26个数据文件组成,11.0 MB(压缩为1个文件,168 KB),反映2013‒2018年海南岛人口时空分布变化特征。
关键词:海南岛;人口;NPP/VIIRS夜间灯光数据;回归模型;人口空间化
DOI: https://doi.org/10.3974/geodp.2022.01.08
CSTR: https://cstr.escience.org.cn/CSTR:20146.14.2022.01.08
数据可用性声明:本文关联实体数据集已在《全球变化数据仓储电子杂志(中英文)》出版,可获取:
https://doi.org/10.3974/geodb.2021.08.02.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2021.08.02.V1.
夜光遥感影像能够为国民生产总值、人口、电力消费、温室气体排放、贫困指数、基尼系数等社会经济参数的估算提供重要依据[1]。目前普遍使用的夜间灯光数据为DMSP/ OLS夜间灯光数据和NPP/VIIRS夜间灯光数据,均来源于美国国家海洋与大气管理局,两者的区别在于DMSP/OLS夜间灯光遥感影像的时间序列为1992‒2013年,空间分辨率为1 km×1 km,而NPP/VIIRS夜间灯光遥感影像的时间序列为2012年至今,空间分辨率为500 m。同时,NPP/VIIRS夜间灯光数据由Suomi-NPP卫星利用VIIRS拍摄,采用极地轨道,由多幅无云影像拼接得到[2],NPP‒VIIRS传感器有22个波段,波长范围在0.4‒12 μm,涵盖可见光和红外光谱区间,光谱分辨率为16 bit,空间分辨率为500 m,更宽波段的辐射探测计和在轨辐射校正技术有效地提高了灯光影像的质量[3]。NPP/VIIRS夜间灯光数据得益于其更强的灯光捕获敏感度、更高的空间分辨率与时间分辨率,更适合中小范围的人类活动研究[4],基于此提取的人口空间分布信息能够直观地反映人口的时空演变特征。因此本数据集基于NPP/VIIRS夜间灯光数据进行处理,并与海南岛18个市县的人口统计数据进行相关性分析回归建模得到海南岛人口空间分布信息,以期为海南省国土空间规划或科学研究提供数据基础及决策依据。
《海南岛夜间灯光与统计数据融合的人口500-m栅格数据集(2013‒2018)》[5]的名称、作者、地理区域、数据年代、时间分辨率、空间分辨率、数据集组成、数据出版与共享服务平台、数据共享政策等信息见表1。
| 条目 | 描述 |
|---|---|
| 数据集名称 | 海南岛夜间灯光与统计数据融合的人口500-m栅格数据集(2013‒2018) |
| 数据集短名 | PopulationHaiNan_2013-2018 |
| 作者信息 | 赖先龙,海南大学,1309842628@qq.com 韩念龙,海南大学,nlhan@hainanu.edu.cn 黄鹏锦,海南大学,1538541468@qq.com |
| 地理区域 | 海南岛 |
| 数据年代 | 2013‒2018年 |
| 时间分辨率 | 逐年 |
| 空间分辨率 | 500 m |
| 数据格式 | .tif |
| 数据量 | 11.0 MB(压缩后168 KB) |
| 数据集组成 | 本数据集由海南岛2013‒2018年夜间灯光与统计数据融合的人口500-m栅格数据集和海南岛2013‒2018年间增长人口空间分布500-m栅格数据组成 |
| 基金项目 | 海南省社科重大项目(HNSK(ZD)19-119));海南省自然科学基金(2019RCO16);海南大学(kyqdsx1962) |
| 数据计算环境 | ArcGIS |
| 出版与共享服务平台 | 全球变化科学研究数据出版系统 http://www.geodoi.ac.cn |
| 地址 | 北京市朝阳区大屯路甲11号100101,中国科学院地理科学与资源研究所 |
| 数据共享政策 | (1)"数据"以最便利的方式通过互联网系统免费向全社会开放,用户免费浏览、免费下载;(2)最终用户使用"数据"需要按照引用格式在参考文献或适当的位置标注数据来源;(3)增值服务用户或以任何形式散发和传播(包括通过计算机服务器)"数据"的用户需要与《全球变化数据学报(中英文)》编辑部签署书面协议,获得许可;(4)摘取"数据"中的部分记录创作新数据的作者需要遵循10%引用原则,即从本数据集中摘取的数据记录少于新数据集总记录量的10%,同时需要对摘取的数据记录标注数据来源[6] |
| 数据和论文检索系统 | DOI, CSTR, Crossref, DCI, CSCD, CNKI, SciEngine, WDS/ISC, GEOSS |
本文研究区域为海南岛陆域范围(不含三沙市)的18个市县,分别是海口市、三亚市和儋州市3个地级市以及15个省直辖县。海南岛位于东经108°37′‒111°03′,北纬18°10′‒20°10′之间,全岛地势四周低平,中间高耸,呈穹隆山地形。
本数据集的原始夜间灯光数据主要是2013‒2018年间的NPP/VIIRS夜间灯光遥感月度合成均值影像,来源于美国国家海洋和大气管理局(NOAA)中的NPP/VIIRS传感器数据集[1],分辨率为500 m;海南岛行政区划数据来源于中国自然资源数据库[2];海南岛的人口统计数据来自海南省统计年鉴[7]。
本研究技术路线包括NPP/VIIRS夜间灯光数据获取,数据处理以及人口空间化建模三部分(图1)。首先是原始数据的获取;其次是数据处理,虽然NPP/VIIRS夜间灯光数据的月度数据剔除了闪电、月光照射及云覆盖的影响,但仍然存在极光、火光、船只灯光等短暂性光源和背景噪声的影响,因此需要将原始数据进行处理得到稳定的年度均值数据,因此NPP/VIIRS夜间灯光数据处理包括年度均值影像合成、消除负值和极高值、去除不稳定光源以及连续校正等过程[8,9]。
图1 数据研发技术路线图
下载人口空间化建模包括相关性分析、回归建模及修正等内容。首先以市县为单位,对2013‒2017年的各市县夜间灯光亮度值与相对应的人口统计数据进行相关性分析。进而对人口统计数据和夜间灯光亮度值进行回归建模分析。结果显示三次线性函数模型的R2最高,因此采用三次线性函数模型。利用2018年数据验证上述模型的人口模拟精度,结果表明,海南岛整体模拟人口平均误差为36.55%,并存在个别市县误差远大于平均误差。为进一步提高人口模拟的精度,将海南岛18个市县按照误差大小进行分类分别二次建模并检验模型精度,结果表明,总体误差为18.19%,人口模拟精度得到了显著提高。最后,由于模拟的人口数量存在一定的误差,因此引入修正系数要进一步消除误差,得到最靠近真实的人口空间分布信息。
2013‒2018年海南岛人口空间化数据集由26个数据文件组成,数据命名采用"主题+时间"的形式命名,如P2013。
将原始数据经过上述技术路线处理后得到分辨率为500 m×500 m的《海南岛夜间灯光与统计数据融合的人口500-m栅格数据集(2013‒2018)》,其中数据单位为500m栅格上的人口总量。由数据集得到的2013‒2018年海南岛人口数量时空分布规律表明,海南岛人口分布整体上较为分散,城区人口密集,而乡镇人口稀疏,沿海市县人口多,而中部市县人口少,形成北海口、南三亚的两极分布特征(图2、图3)。在2013‒2018年间增长的人口主要分布在各市县的城区外围,而城区中心人口增长较少,如海口、三亚两市该情况尤为明显,整体上增长人口主要分布在沿海市县(图4)。
图2 2013年海南岛人口空间化
下载
图3 2018年海南岛人口空间化
下载
图4 2013‒2018年人口增长
下载将2013‒2018年的NPP/VIIRS夜间灯光数据通过模型反演得出历年人口数量,并依据统计数据验证反演的误差。经过误差分类进行的二次建模人口模拟精度较第一次建模人口模拟精度结果有显著提高,回归系数R2均为0.9以上,与曹丽琴等基于线性模型估算湖北地级市人口[11](R2为0.923)以及李翔基于空间回归模型估算上海市各乡镇街道的人口[10](R2为0.91)较为接近。模型反演结果的平均相对误差为18.19%,精度与赵宇等基于NPP/VIIRS夜间灯光数据构建线性回归模型估算大连市各辖区人口的平均相对误差为18%相近[12],高于王璇璇等基于NPP/VIIRS夜间灯光数据构建空间滞后模型估算江苏省各县市常住人口,平均相对误差为26.09%[13]和陈晴等基于DMSP/OLS夜间灯光数据构建线性回归模型估算黄河三角洲高效生态经济区19个县市区人口,平均相对误差为30%[14]。
在市县层面,各市县的反演结果中五指山市的误差较大,原因在于五指山市位于海南岛中部山区腹地,城市地形起伏较大且人口较为稀疏,单凭夜间灯光数据难以准确反演其人口空间分布数量。根据海南省行政区划,各市县按照地级市、县级市、普通县和少数民族自治县划分,可以发现二次建模的结果中地级市的人口模拟误差最小,其次是少数民族自治县,县级市误差略高于少数民族自治县,一般普通县的误差最大,不同行政规模的市县之间的人口模拟精度并没有规律性区别。
同时,付晶莹等人针对我国人口空间分布特征及区域差异进行人口区划,筛选出与人口分布相关性强的土地利用类型作为建模参数,建立多元回归模型实现人口空间化,并利用城市人口密度、交通状况、DEM及总量控制四个因子对模型进行校正[15],最终得到全国人口分布数据集研究产品,其中海南岛2010年人口空间化数据的平均相对误差仅为0.92%[16],因此,未来可参考这些方法结合夜间灯光数据进行人口空间化以提升精度。
| 类别 | 市县 | 统计人口(万人) | 一次模拟人口(万人) | 误差 | 二次模拟人口(万人) | 误差 |
|---|---|---|---|---|---|---|
| 地级市 | 海口市 | 230.23 | 258.12 | 12.12% | 194.61 | 15.47% |
| 三亚市 | 77.39 | 121.09 | 56.46% | 74.51 | 3.72% | |
| 儋州市 | 99.84 | 90.34 | 9.52% | 117.89 | 18.08% | |
| 平均相对误差 | 26.03% | 12.42% | ||||
| 县级市 | 文昌市 | 56.89 | 42.34 | 25.58% | 51.57 | 9.35% |
| 琼海市 | 51.57 | 43.34 | 15.97% | 52.66 | 2.11% | |
| 万宁市 | 57.86 | 47.28 | 18.28% | 57.13 | 1.27% | |
| 五指山市 | 10.71 | 26.27 | 145.25% | 16.36 | 52.77% | |
| 东方市 | 42.97 | 62.63 | 45.74% | 51.95 | 20.89% | |
| 平均相对误差 | 50.16% | 17.28% | ||||
| 普通县 | 定安县 | 29.76 | 30.05 | 0.99% | 39.55 | 32.88% |
| 屯昌县 | 26.85 | 24.67 | 8.11% | 35.09 | 30.68% | |
| 澄迈县 | 49.44 | 57.02 | 15.33% | 69.19 | 39.94% | |
| 临高县 | 45.1 | 35.25 | 21.84% | 44.32 | 1.73% | |
| 平均相对误差 | 11.57% | 26.31% | ||||
| 少数民族自治县 | 乐东县 | 48.27 | 48.94 | 1.39% | 59.08 | 22.39% |
| 琼中县 | 18.02 | 24.86 | 37.95% | 14.87 | 17.49% | |
| 保亭县 | 15.28 | 26.68 | 74.59% | 16.8 | 9.93% | |
| 陵水县 | 33.39 | 52.65 | 57.67% | 43.06 | 28.95% | |
| 白沙县 | 17.34 | 26.71 | 54.05% | 16.83 | 2.91% | |
| 昌江县 | 23.35 | 36.68 | 57.07% | 27.26 | 16.75% | |
| 平均相对误差 | 47.12% | 16.40% | ||||
| 整体平均相对误差 | 36.55% | 18.19% | ||||
NPP/VIIRS夜间灯光数据人口空间化建模能基本上反映人口的空间分布,但是五指山市的模拟精度较差。其原因可能是由于五指山市位于中部山区,地形起伏度较大,使得夜间灯光数据受较大影响,且该市的人口数量与其他市县差距较大[14]。针对该问题,未来将考虑对研究区进行更详细的人口区划,同时结合土地利用数据,POI等多源数据[17]并采用多地理因子加权法、空间加权回归或神经网络等模型以提高反演精度。
本数据集是基于NPP/VIIRS夜间灯光数据反演的500-m分辨率海南岛人口空间化数据,该数据集能够反映出海南岛2013‒2018年的人口时空动态分布情况,可为当前的海南国土空间规划提供辅助依据。
作者分工:韩念龙对数据集的开发做了总体设计;赖先龙采集和处理了夜间灯光数据数据;赖先龙设计了模型和算法;黄鹏锦做了数据验证;赖先龙撰写了数据论文等。
利益冲突声明:本研究不存在研究者以及与公开研究成果有关的利益冲突。
[1] 美国国家海洋和大气管理局. https://ngdc.noaa.gov.
[2] 中国自然资源数据库. http://www.naturalresources.csdb.cn/index.asp.